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Abstract 

The problem is addressed of denning the values of functions, whose variables tend to 
infinity, from the knowledge of these functions at asymptotically small variables close 
to zero. For this purpose, the extrapolation by means of different types of self-similar 
approximants is employed. Two new variants of such an extrapolation are suggested. 
The methods are illustrated by several examples of systems typical of chemical physics, 
statistical physics, and quantum physics. The developed methods make it possible to 
find good approximations for the strong-coupling limits from the knowledge of the 
weak-coupling expansions. 
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1 Introduction 



In a variety of problems of chemistry and physics [1], there is the necessity of defining 
the properties of systems in the strong-coupling limit, when, however, this limit cannot 
be accessed directly, but when, because of the complexity of the problem, only the weak- 
coupling expansion is available. Then the principal question is whether one can infer the 
strong-coupling properties from a weak-coupling expansion? 

In the present paper, we address this principal question and show that the problem can 
be resolved by employing the self-similar approximation theory [2-12] allowing for the extrap- 
olation of asymptotic series to finite values of expansion parameters. Here we concentrate on 
the question of how to better organize an extrapolation not merely to finite parameter values, 
but to their infinite values. We show that the self-similar extrapolation is a convenient tool 
for this purpose. Since this approach allows for the construction of different approximation 
schemes, it is useful to compare the latter in order to decide which of them better suits the 
posed problem. It is also necessary to make some modifications for adjusting the methods 
to the case, when the behavior at infinity is of interest. Here we make such modifications 
and suggest two new ways of constructing self-similar approximants, which have not been 
analyzed in our previous publications. 

To illustrate the applicability of the developed methods and to compare their accuracy, 
we consider several problems whose mathematical structure is typical of chemical physics, 
quantum physics, and statistical physics. We compare different variants of the self-similar 
approximants. Also, we consider the Pade approximants [13] and show that the accuracy of 
the best Pade approximant is always inferior to that of the best self-similar approximant. 
In some cases, the Pade approximants are not applicable at all, yielding qualitatively wrong 
results, while there always exists an accurate self-similar approximant. 

In Sec. 2, we give the general formulation of the problem, defining the strong-coupling 
exponents and amplitudes. We introduce two novel methods of constructing the self-similar 
approximants yielding the corrected self-similar approximants and the iterated root approxi- 
mants. For the beginning, we choose, in Sec. 3, the Debye-Huckel function from the theory of 
strong electrolytes. The quartic anharmonic oscillator of Sec. 4 provides the weak-coupling 
asymptotic series with the mathematical structure that is common for many problems of 
quantum mechanics and field theory. In Sec. 5, we consider the expansion factor of a poly- 
mer chain. The ground-state energy for a model of one-dimensional string is calculated in 
Sec. 6. A more realistic model of fluctuating fluid membrane is treated in Sec. 7. A dilute 
gas of strongly interacting fermions is studied in Sec. 8. Then, in Sec. 9, the strong-coupling 
limit of a one-dimensional Bose system is analyzed. The problem of calculating the Bose- 
Einstein condensation temperature of an interacting Bose gas is investigated in Sec. 10. 
Finally, Sec. 11 concludes. 

Some of the examples, analized in the present paper, have been considered earlier [14,15], 
however, by other methods. The main difference of this paper from our previous publications 
is threefold: (i) We consider a large variety of problems to show the generality of our theory, 
(ii) We give a comparative analysis of different approximants. (iii) We suggest new variants 
of self-similar approximants, which have not been considered previously. 
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2 Formulation of general problem 



Let us be interested in a real function f(g) of a real coupling parameter g, when g — > 
oo, which corresponds to the strong-coupling limit. But assume that the problem is so 
complicated that we are not able to calculate this limit directly. Instead, we can get only 
the weak-coupling expansion, when perturbation theory is applicable. Then, we can find the 
value 

/(<?) * fk(g) (g o) (l) 

for asymptotically small g, which is given by the weak-coupling expansion 

k 

f k (g) = fo(g)Y. a -9 n ■ ( 2 ) 

n=0 

Here the factor fo(g) is not expandable in powers of g. Without the loss of generality, we 
can always set a — 1, which, however, is not principal. 



2.1 Strong-coupling exponents and amplitudes 

^From some additional information, we may know that the behavior of the sought function 
at large coupling parameter is as follows: 

/(<?) ~ Bgf* (g^oo). (3) 

The power index j3 is called the strong- coupling exponent and can be defined as 

^hm^M. (4) 
g^oo In g 

This can also be written as 

_ .. dlnf(g) g df(g) 

f3 = km JKyJ = km -f- -^LL . 
g^oo ding g^oo f(g) dg 

We assume tkat tkis exponent is known, and wkat we need to find is tke strong- coupling 
amplitude 

B = lim (5) 

g^oo gP 

Tke problem we confront is kow to find tke strong-coupling amplitude (5) from tke knowledge 
of tke weak-coupling expansion (2)? 



2.2 Corrected self-similar approximants 

An effective extrapolation can be done by means of tke self-similar approximation tkeory [2- 
12]. Tken, from tke weak-coupling expansion fk(g), we construct a self-similar approximant 
fk(g)- Wken tke bekavior of tke sougkt function is supposed to be exponential, we kave to 
use tke self-similar exponential approximants [16,17]. But in tke considered case (3) of tke 
power law bekavior, it is better to employ tke self-similar factor approximants [18-21]. 
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Suppose, we have found an approximant f k {g). This can be an exponential approximant 
[16,17] or a factor approximant [18-21]. Now we suggest a novel way of correcting this 
approximant and, respectively, improving the accuracy of calculating the strong-coupling 
amplitude. 

Let an approximant f k (g) be given, whose parameters are defined by the accuracy- 
through-order procedure, which implies that the approximant f k (g) is re-expanded in powers 
of g and compared, term by term, with the initial expansion, so that 

/*(<?) ^ /*(<?) (g->0). (6) 

This is also called the re- expansion procedure. More details, describing this procedure, can 
be found in Ref. [15]. 

In constructing the approximant f k (g), we keep in mind the limiting condition 

lim MM = /3. (7) 
g-^oo In g 

Our aim is to find the strong-coupling amplitude 

B k = lim Mf. (8) 



Now, let us introduce a correcting function 



C k+P (g) - ^) (g^ 0) , (9) 

which is defined as an expansion in powers of g, such that 

k+p 

C k+p (g) = J2 b "9 n i (10) 

n=0 

whose coefficients b n are, evidently, functions of the initial coefficients a n from series (2). 
Since the parameters of f k {g), according to Eq. (6), are defined through the re-expansion 
procedure, the correcting function (10) contains not all powers of g but only those starting 
from the order k + 1, that is, 

k+p 

C k+P (g) = l+ b n9 n - (11) 

n=k+l 

This fact essentially simplifies the construction of a self-similar approximant C% + (g) from 
the correcting function (11). Constructing the approximant C% + (g), we impose the limiting 
condition 

lim C* k+ Jg) = const . (12) 

g^oo 

The corrected self-similar approximant is defined as 

lk +P {g) = r k (g)C* k+p (g) . (13) 
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Because of condition (12), 



9^00 In g 

Therefore the strong-coupling exponent does not change, 



lim ^±M = . (14) 



lim '-{^ M = lim ^M = fi. (15) 



g^oo In g g^oo In g 

But the corrected strong- coupling amplitude 



changes to 



where B k is given by Eq. (8). 



B k+P ee lim f -^± (16) 

9->oo gP 



B k+P = B k lim C* k+p (g) , (17) 



2.3 Iterated root approximants 

In addition to self-similar exponential approximants [16,17] and factor approximants [18- 
21], there are self-similar root approximants [10-12]. The latter serve as a convenient tool 
for interpolating functions whose strong-coupling behavior is described by means of known 
strong-coupling expansions, when the root approximants can be uniquely determined [22]. 
However, the parameters of the root approximants cannot be uniquely defined through the 
re-expansion procedure in the region of weak coupling. Here we suggest a method allowing 
us to derive the parameters of the root approximants by means of the re-expansion procedure 

(6) . ' 

The idea of the method, allowing for a unique derivation of the parameters of the root 
approximants, is as follows. In order to remove the multiplicity of solutions, arising in the re- 
expansion procedure applied to the general form of the root approximants, it is necessary to 
impose some restrictions on the definition of the related parameters. Thus, we can keep the 
same parameters of the lower-order approximants in the higher-order approximants, leaving 
there unknown only the highest-order parameter that is to be found from the re-expansion 
procedure. And the highest-order power has to be such that to satisfy the limiting condition 

(7) . Since in this method, the lower-order root approximants are inserted into the higher- 
order approximants, the resulting expressions can be called, iterated root approximants. 

To be concrete, let us assume that the factor f (g) possesses the strong-coupling behavior 

as 

fo(g) ~ Ag a (g^oc). (18) 
The first-order root approximant 

R 1 (g) = fo(g)(l + A 1 g)\ (19) 

where we set 

7 = P ~ a , (20) 
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is the same as the first-order factor approximant and it is uniquely defined with A\ = ai/7. 
In the second-order root approximant 

Riig) = foig) ((i + A l9 f + A 29 2 y' 2 , (21) 

we keep the same A\ as before, and 7 is given by Eq. (20). Hence, we need to find only A2 
from the re-expansion procedure. In the third-order root approximant 

Rs(g) = fo(9) (((1 + A l9 f + A 2 g 2 f 2 + A 3 g 3 Y /3 , (22) 

the parameters Ai and A 2 are kept the same as in the previous expression (21), while A 3 is 
defined through the re-expansion procedure. Such an iteration construction continues, with 
the general /c-approximant being 

R k (g) = foig) ((• • • (i + Mgf + A 2 g 2 f 2 + ... + A k9 k y /k , (23) 

where all parameters A n , with n — 1,2, . . . ,k — 1 are the same as in the k — 1-order approx- 
imant and the parameter A k is defined by the re-expansion procedure. 

Constructing in this way the iterated root approximants yields the strong-coupling am- 
plitudes 

B k =lim^M. (24) 

In particular, 

B 1 — AAJ , B 2 = A (A 2 , + A 2 ) l/2 , B 3 = A [{A\ + A 2 ) 3/2 + A^^ , 
and so on, with the A;-order amplitude 

B k = A^...(Al + A 2 f 2 + A 3 y /3 + ... + A k y / . (25) 

An accurate evaluation of the strong-coupling amplitude is the main aim of the suggested 
scheme. In the following sections, we shall consider several examples of explicit calculations 
of strong-coupling amplitudes, comparing the accuracy of different approaches. 



3 Debye-Hiickel function 

Let us start with the case, when the sought function is actually known. This will allow us 
to easily determine the accuracy of each scheme. Let us consider the Debye-Hiickel function 
that is met in the theory of strong electrolites [23,24]. In dimensionless units, the function 
reads as 

M = -g ~ J2 (1 " • (26) 

In view of Eqs. (2)-(5), we have fo(g) — 1, /3 — —1, and 

B = 2 . (27) 
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The coefficients, corresponding to the weak-coupling expansion (2), are 



ai ~ ~ 3 ' ° 2 " 12 ' ° 3 " " 60 ' 
111 1 



360 2520 20160 181440 

The Pade approximants do not provide good accuracy. The best of them, 

p / \ _ 1 + gig + C 2 gf 2 + C 3 q 3 
J/4WJ " 1 + C 4 <? + C 5 <? 2 + C 6 <f + > 

possessing the correct strong-coupling exponent /3 = — 1, and invoking all seven expansion 
coefficients a n , gives the strong-coupling amplitude 

£3/4 = lim gPy 4 (g) = — = 1-6 , 

g-S-oo C7 

with an error e(5 3 / 4 ) = —20%. 



3.1 Self-similar factor approximants 

The factor approximant, using the expansion terms up to the third order, that is, involving 
the coefficients 01,02,03, reads as 

K(g) = (l + A 1 g) ni (l + A 2 g)'» , (28) 

with the condition 

rai + n 2 = -l. (29) 
This gives the strong-coupling amplitude 

B 3 = A 7 ? A™ 2 = 1.64 , (30) 

whose error is e(B 3 ) = —18%. 

The fifth-order factor approximant 

/*(<?) = (1 + A l9 r(i + A 2g ra + A 39 r , (31) 

invoking the weak-coupling expansion terms up to the fifth order in g and satisfying the 
condition 

n\ + n 2 + n 3 — -1 , (32) 
yields the strong-coupling amplitude 

B 5 = A^A^Af = 2.295 , (33) 

with an error e(B 5 ) = 15%. 

Note that the parameters of each of the factor approximants are defined through the 
re-expansion procedure. So, the values of and in each order are different. We use the 
same letters in different approximants just for the simplicity of notation. 
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The factor approximant, using all seven expansion terms, is 

r 7 { 9 ) = (i + A l9 r(i + A 29 r(i + ^ra + A 49 r , (34) 

under the condition 

n\ + n 2 + n 3 + n 4 = -1 . (35) 
This gives the strong-coupling amplitude 

B 7 = A^A^A^A^ 4 = 1.799 , (36) 

with an error e(B-j) = —10%. 

3.2 Corrected factor approximants 

Following the scheme of Sec. 2, we can define the corrected factor approximants. To this 
end, let us introduce the correction function 

which yields 

C 7 (g) = 1 + hg 4 + b 5 g 5 + b 6 g 6 + b 7 g 7 . (38) 
^From this expansion, we construct the factor approximant 

C* 7 (g) = 1 + b 4 g 4 (l + D 39 r*(l + D 4 g) ni , (39) 

under the condition 

4 + n 3 + n 4 = . (40) 
The corrected factor approximant is defined as 

frig) = mg)c* 7 (g) . (41) 

This yields the strong-coupling amplitude 

B 7 = B 3 (l + b 4 D^D2 4 ) = 1.944 , (42) 
with an error e(B 7 ) = —2.8%. 

3.3 Simple root approximants 

A simple root approximant 

r* 3 (g) = ((l + A l9 r+A 2 g 2 )- 1/2 , (43) 

in which the parameters are defined trough the re-expansion procedure, gives the strong- 
coupling amplitude 

B 3 = -^= = 2.164, (44) 
v Ao 



whose error is e(B 3 ) = 8.2%. 



3.4 Corrected root approximants 

The above simple root approximant can be corrected according to Sec. 2. For this purpose, 
we define the correction function 

Cr(9) - ^ (9^ 0) , (45) 

corresponding to the expansion in powers of g up to the seventh order, and construct the 
root approximant 

C* 7 (g) = 1 + M 4 ((1 + D 39 r* + ZVr 2 • (46) 
Then the corrected root approximant is 

r*r(9) = r* 3 (g)C* 7 (g) . (47) 

The latter gives the strong-coupling amplitude 



1 (. h 



B 7 = M + *) = 1.875, (48) 



with an error e(B 7 ) = —6.3%. 

3.5 Iterated root approximants 

The iterated root approximants, according to Sec. 2, have the form 



Rk 



k (g) =((... (1 + A l9 f + A 2 g 2 f 2 + ... + A k g k ) m , (49) 



where (3 — — 1. Therefore, the strong-coupling amplitude of /c-order is 



B k =((...(Al + A 2 f 2 + A 3 y /3 + ... + A k )j' 



(50) 



/,From here, we find 

B 1= 3, .02 = 2.449, B 3 = 2.229, B 4 = 2.127 , 

5 5 = 2.067 , B 6 = 2.032 , B 7 = 2.009 . 

The seventh order amplitude has an error e(B 7 ) = 0.45%. Thus, the iterated root approx- 
imant (49) of the seven-th order is the most accurate among all above considered approxi- 
mants. 

4 Quartic anharmonic oscillator 

The anharmonic oscillator, with the Hamiltonian 



is the model imitating many systems in quantum chemistry, atomic physics, condensed- 
matter physics, and field theory. Here, g is the coupling parameter and we set x G (—00, 00) 
and g G [0, 00). The ground-state energy for this Hamiltonian, found by means of the 
weak-coupling perturbation theory [25,26], when g — > 0, reads as 

e(g) ~ a + a x g + a 2 g 2 + a 3 g 3 + a A g 4 + a 5 g 5 + a 6 g 6 + a 7 g 7 , (52) 

with the coefficients 

1 3 21 333 

a — 77 ? a l — T ? «2 — — j O3 



2 ' 1 4 ' z 8 ' 3 16 ' 
30885 916731 65518401 2723294673 

a 4 — 7^7T~ ) a 5 — — ^7 — j OjQ — — — , a-j 



128 ' 256 ' 1024 ' 2048 

The strong-coupling limit for the ground-state energy yields 

e(g) ~ 0.667986 g 1/3 (g ->■ 00) . (53) 

Hence, the strong-coupling exponent /3 = 1/3 and the strong-coupling amplitude is 

B = 0.667986 . (54) 

As is seen, the coefficients of the weak-coupling expansion quickly grow, and perturbation 
theory is divergent for any nonzero g. The standard Pade approximants for this problem are 
not applicable at all, being unable to satisfy the fractional strong-coupling exponent j3 = 1/3. 
It is possible to introduce the power-transformed approximants, similarly to Ref. [14], fixing 
the exponent /3 = 1/3. But this way does not provide a great accuracy improvement. The 
best power-transformed Pade approximant gives an error of 8%. 

4.1 Self-similar factor approximants 

Using the third order of the weak-coupling expansion (52), we get the self-similar factor 
approximant 

e;(g) = ^l + A ig r(l + A 2 gr , (55) 

under the condition 

1 . . 

ni + n 2 = - . (56) 

This gives the strong-coupling amplitude 

B 3 = X - AT AT = 0.75 , (57) 

with an error e(B 3 ) = 12%. 

The fifth-order factor approximant has the form 

et(g) = 1(1 + A 1 g) ni (l + A 29 r(l + Asg)** , (58) 

with the condition 

1 . . 

n\ + n 2 + n 3 = - . (59) 
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Then the strong-coupling amplitude is 



5 5 = i A^AfAf = 0.725 , 

whose error is e(B 5 ) = 8.5%. 

The seventh-order factor approximant is 

e* 7 (g) = i(l + A l9 r(l + A 29 r(l + A 3 g)^{\ + A i9 f 
with the condition 

1 

ni + n 2 + n 3 + n A — - . 
^From here, we have the strong-coupling amplitude 

1 



B 7 = - A?A?AZ 3 Ap = 0.712 , 



with an error e(B 7 ) = 6.6%. 



4.2 Corrected factor approximants 

Defining the correction function 

we get the series 

7 

C 7 (g) = l + J2 b ^ n ■ 

n=4 

The latter generates the factor approximant 

C* 7 (g) = 1 + Kg\\ + D 3 g) n \\ + D A g) ni , 

where 

4 + n 3 + n 4 = . 

The corrected factor approximant 

e 7 {g) = e* 3 (g)C* 7 (g) 

yields the strong-coupling amplitude 

B 7 = I A^Af (1 + b 4 D^ 3 D2 4 ) = 0.728 , 
with an error £(£7) = 9%). 
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4.3 Simple root approximants 

Simple root approximants do not give good accuracy. For instance, the root approximant 

rt(g) = \ (((1 + A l9 r + A 2 g 2 ) 2 + A^) 1 " (70) 
results if the strong-coupling amplitude 

B 5 = \ Af = 0.824 , (71) 
whose accuracy is characterized by the error e(B 5 ) = 23%. 

4.4 Iterated root approximants 

According to Sec. 2, the iterated root approximants are given by 

Ri(g) = \{l + A l9 )^ , R 2 {g) = \{{l + A l9 f + A 2 g 2 ) 1/6 , 

R3(g)= 1 -(((l + A l9 f + A 29 *) 3/2 + Al) 1/9 , (72) 

and so on. However, the approximant R^(g) is complex, and we limit ourselves by Eqs. (72). 
The related strong-coupling amplitudes are 

B 1 = \ A{ /3 = 0.825 , B 2 = \(A\ + A 2 ) 1/6 = 0.572 , 

Bz = \({A\ + A 2 f 2 + A^ 9 = 0.855 . (73) 
Their errors are e{B 1 ) = 24%, e(B 2 ) = -14%, s(B 3 ) = 28%. 

4.5 Power-transformed approximants 

It is possible to define the power transformation [14] as 

P k (g,m)^f™(g) (g -+ 0) , (74) 

producing the series 

k 

P k (g,m) = J2 b n(m)g n ■ 

n=0 

Then, constructing a self-similar approximation on the basis of this series, we get P k (g,m). 
The final answer is given by the inverse power transformation 

fk(9) = [Pk(9,™)] 1,m ■ (75) 

This way was analized in Ref. [14], where the transformation power m was defined by 
variational optimization. Such a way, however, is rather cumbersome and yields the accuracy 
improvement not better than the simpler methods considered here. 

Among the methods, studied in this Section, the self-similar factor approximants provide 
the best accuracy for the strong-coupling amplitude. 
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5 Polymer expansion factor 

The expansion factor a(g), as a function of the coupling parameter g, for a polymer chain 
with excluded interactions, can be calculated by means of perturbation theory [27] resulting 
in the series 

a(g) = 1 + aig + a 2 g 2 + a 3 g 3 + a 4 g A + +a 5 g 5 + a 6 g 6 , (76) 
with the coefficients 

ai = - , a 2 = -2.075385396 , a 3 = 6.296879676 , 

a 4 = -25.05725072 , a 5 = 116.134785 , a 6 = -594.71663 . 
Numerical fitting [27,28] gives the phenomenological formula 

a(^)= (1 + 7.524^ + 11.06/) °- 1772 . (77) 

This implies that the strong-coupling exponent (3 = 0.3544, which is in agreement with other 
numerical simulations [29], where the exponent 

was calculated, giving v = 0.5877. Therefore the strong-coupling amplitude is 

B = 1.531 . (79) 

5.1 Self-similar factor approximants 

The self-similar factor approximants start with the lowest order 

a* 1 (g) = (l + A 1 g f , (80) 

in which A\ = a\j (3 and j3 = 0.3544. Then the strong-coupling amplitude is 

B 1 = A{ = 1.599 , (81) 

with an error e(Bi) = 4.4%. Increasing the approximation order improves the accuracy. 
Thus, the fifth-order approximant 

a* 5 (g) = (1 + A l9 r(l + A 29 r(l + A.g)^ , (82) 

where 

wi + n 2 + n 3 = /3 = 0.3544 , (83) 
gives the strong-coupling amplitude 

B 5 = A\ x A\ 2 Al z = 1.541 , (84) 

with an error e(B 5 ) = 0.65%. 
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5.2 Corrected factor approximants 

Introducing the correction function 

CM * ^ (</ 0) , (85) 

and constructing the related factor approximant 

C;(g) = l + b 2 g 2 (l + D ig )- 2 , (86) 
we come to the corrected factor approximant 

a 3 (g) = al(g)C;(g) . (87) 
This gives the strong-coupling amplitude 

B 3 =A^(l + ^j= 1.552 , (88) 

with an error e(B 3 ) = 1.4%. 

The accuracy of the corrected factor approximants is close to that of the directly con- 
structed self-similar factor approximants. 

5.3 Simple root approximants 

Simple root approximants provide sufficiently good accuracy. Thus the approximant 

r*M = (((1 + A l9 r + A 2 gT + A^) m (89) 
yields the strong-coupling amplitude 

B 5 = A% /3 = 1.537 , (90) 

with an error e(B 5 ) = 0.4%. 

5.4 Corrected root approximants 

For the correction function 

C,(g)^^f\ (<7->0), (91) 
the corresponding self-similar root approximant reads as 

C*(g) = 1 + b 2 g 2 ((1 + Dig) 1711 + D 2 g 2 ) 1 . (92) 
The corrected root approximant 

r 5 (g) = al(g)C;(g) (93) 
results in the strong-coupling amplitude 

B 5 = Ai(l + ^-) = 1.542, (94) 



whose error is e(B 5 ) = 0.72%. 

The simple root approximants, being already sufficiently accurate, are not improved by 
the corrected variants. 
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5.5 Iterated root approximants 



For the iterated root approximants, following the scheme of Sec. 2, we find the strong- 
coupling amplitudes 

B 1 = 1.599 , B 2 = 1.544 , B 3 = 1.549 , 

B A = 1.539 , B 5 = 1.541 , B Q = 1.537 . (95) 

The error of the highest-order approximant is e(B 6 ) = 0.4%. 

All variants of the approximants, studied in the present Section, are of about the same 
accuracy. We have also considered the power-transformed approximants [14] , whose accuracy 
is found to be close to that of the iterated root approximants. 



6 Oscillating fluid string 

In biological and chemical applications, there exists an important class of membranes called 
fluid membranes [30]. The model of a one-dimensional fluid string, oscillating between two 
walls, was advanced by Edwards [31]. It has been shown [31-33] that the free energy of the 
string coincides with the ground-state energy of a quantum particle in a one-dimensional 
box. This energy, as a function of a finite wall stiffness g, reads as 

^) = ^|l + ^ + 7\/l + ^l • 06) 



8g 2 




In the small- stiffness expansion 



7T 2 



( 97 ) 



n=0 



the coefficients are 



ao = l, a 1 = i, a 2 = l, «3 = ^, 

1 1 
a4 = , 05 = — -77-7— , a6 = , 07 



131072 ' 16777216 
The rigid walls correspond to the stiffness g — > 00, when 



TT 2 



E(po) = — = B = 0.077106 . (98) 
128 

Hence, the strong-stiffness exponent is (3 = 0. 

Our aim is to approximate the rigid-wall limit (98) employing the small-stiffness expan- 
sion (97). 

Pade approximants are not applicable for this problem, giving negative values of the 
strong-stiffness amplitude. So, we shall apply the self-similar approximants. 
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6.1 Simple root approximants 

The root approximants of low orders are not sufficiently accurate. Thus, the root approxi- 
mant 

r* 3 (g)=E (g)((l + A ig r+A 2 g 2 ) , (99) 
where E (g) = ir 2 /8g 2 , gives the strong-stiffness amplitude 

7T 2 

B 3 = — A 2 = 0.0544 , (100) 

8 

with an error e(B 3 ) = —29%. And the higher-order root approximants are not convenient, 
being not uniquely defined through the re-expansion procedure. 



6.2 Corrected root approximants 

As is shown in Sec. 2, we can introduce the corrected root approximants. For instance, the 
correction function 

CM - ^ (^o) (ioi) 

generates the self-similar root approximant 

C* 7 (g) = 1 + M 4 ((1 + D ig )^ + D 2 g 2 ) 2 . (102) 
Then the corrected root approximant 

r* 7 (g)=r* 3 (g)C* 7 (g) (103) 

yields the strong-stiffness amplitude 

B 7 = yA 2 (l + |)= 0.0771 , (104) 

with an error e(B 7 ) = —0.02%. 

We may note that the sixth-order corrected root approximant is also of good accuracy, 
with an error B 6 = —0.66%. 



6.3 Iterated root approximants 

The uniquely defined iterated root approximants, introduced in Sec. 2, result in the strong- 
stiffness amplitudes 

B x = 0.019 , B 2 = 0.039 , B 3 = 0.051 , 

B A = 0.058 , B 5 = 0.062 , B G = 0.065 , B 7 = 0.067 . (105) 

The errors of the highest two approximations are e(B 6 ) = —16% and e(B 7 ) = —13%, re- 
spectively. 

Another variant was considered in Ref. [15], invoking the change of the variable g, and 
giving the strong-stiffness amplitude of seventh order with an error of 5%. The best accuracy 
is provided by the corrected root approximant (103) yielding amplitude (104). 
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7 Fluctuating fluid membrane 



For a three-dimensional fluctuating fluid membrane, there are no exact results. Thus, the 
pressure of this membrane is calculated by means of perturbation theory [34] yielding ex- 
pansions 

2 ^ 
y n=0 

The coefficients are known [34] up to the sixth order: 

a = 1 , ai = ^ , a 2 = , a 3 = 2.176347 x 1(T 3 , 

a 4 = 0.552721 x 1(T 4 , a 5 = -0.721482 x 10" 5 , a 6 = -1.777848 x 10~ 6 . 
The strong-stiffness limit has been found by Monte Carlo simulations [35], giving 

p(oo) = 0.0798 ± 0.0003 . (107) 

This implies that the strong-stiffness exponent is (3 = and the strong-stiffness amplitude is 

B = 0.0798 . (108) 

Pade approximants are inapplicable for this problem, resulting in negative values for B. 

7.1 Simple root approximants 

The simple root approximant 

rl{9) = ^ ((1 + A l9 r + A 2 g 2 ) (109) 
gives the strong-stiffness amplitude 

7T 2 

B 3 = — A 2 = 0.056 , (110) 

o 

with an error e{B^) = —30%. The higher-order root approximants are not convenient to 
employ, since they are not uniquely defined. 

7.2 Corrected root approximants 

If we set aj = 0, we can introduce the correction function 

^)-^| (9^0). (Ill) 
The corresponding self-similar root approximant is 

C* 7 (g) = 1 + b A g A ((1 + D 1 g) mi + D 2 g 2 ) 2 . (112) 
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The corrected root approximant is defined as 

prig) = r* 3 (g)C* 7 (g) , (113) 
which gives the strong-stiffness amplitude 



7T 2 , / b. 



!>■; -Mll jjj) = 0.085, (111) 



with an error e{B-j) = 6.5%. 



7.3 Iterated root approximants 

Employing the iterated root approximants of Sec. 2, we have 

B l = 0.019 , B 2 = 0.039 , B 3 = 0.053 , 

B A = 0.061 , B 5 = 0.067 , B Q = 0.071 . (115) 

The error of B 6 is e(B 6 ) = -11%. 

In this way, for this problem, the corrected root approximants are the most accurate 
ones. 

8 Strongly interacting fermions 

The ground-state energy of a diluted Fermi gas, using perturbation theory, can be represented 
[36,37] as an asymptotic expansion 

e(g) ~ a + a x g + a 2 g 2 + a 3 g 3 + a A g 4 (116) 

in powers of the dimensionless coupling parameter 

g = \k F a s \ , (117) 

in which hp is a Fermi wave vector and a s , scattering length. Here g — > and the coefficients 
are 

3 1 

°o = T7T > fl i = _ T~ ' °2 = 0.055661 , 
10 3n 

a 3 = -0.00914 , a 4 = -0.018604 . 
Numerical calculations [38,39] give the strong-coupling limit 

lim e(g) = 0.132 . (118) 

Hence, the strong-coupling exponent is (3 = 0, while the strong-coupling amplitude is 

5 = 0.132. (119) 

Pade approximants are not accurate for this problem, the best of them having an error 
of 30%. 
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8.1 Self-similar factor approximants 

The third-order factor approximant reads as 

e;(g) = ^(l + A ig r(l + A 2 gr , (120) 

under the condition 

ni = n 2 = . (121) 
This yields the strong-coupling amplitude 

B 3 = ^A?A? = 0.174, (122) 

with an error e(B 3 ) = 32%. 

The fourth-order factor approximant can be written as 

e* 4 (g) = 1 + a l9 (l + A^gf^l + A 4 g) n * , (123) 

under the condition 

1 + n 3 + n 4 = . (124) 
Then the strong-coupling amplitude becomes 

Bi = To + aiA z 3A * 4 = °- 162 ' ( 125 ) 

with an error e(B 4 ) = 23%. 

8.2 Corrected factor approximants 

For the correction function 

C*(g) ^ C 4r\ (9^ 0) , (126) 
the factor approximant can be represented as 

Cl{g) = l + b A g\l + D ig )- A . (127) 
The corrected factor approximant 

e 4 (<7) = e* 3 (g)C* 4 (g) (128) 

yields the strong-coupling amplitude 

B A = A A? 1 ^? (i + ^) = °- 1434 > ( 129 ) 
whose error is e(B 4 ) = 8.6%. 
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8.3 Iterated root approximants 

The iterated root approximants, as defined in Sec. 2, exist for the lowest two approximations, 
giving 

Bi = 0.098 , B 2 = 0.169 . (130) 

The higher approximants become complex. The accuracy of B 2 is characterized by an error 
e(B 2 ) = 28%. 

In this way, the best approximation is provided by the corrected factor approximant 
(128), with an error of 8.6%. 

9 One-dimensional Bose system 

The one-dimensional Bose gas with local interactions corresponds to the Lieb-Liniger model 
[40]. The ground-state energy of the gas, in the weak-coupling limit, is described by the 
expansion 

e(g) ~g + a 3 g 3 ^ 2 + a 4 g 2 + a 5 g 5/2 , (131) 

with the coefficients 

a 3 = - — = 0.424413 , a 4 = = 0.065352 , a 5 = -0.017201 . 
In the strong-coupling limit, one gets the Tonks- Girardeau expression 

7T 2 

lim e(g) = — = B = 3.289868 . (132) 

g->oo 3 

More details can be found in the review [41]. 

9.1 Simple root approximants 

The simple root approximant 

r* 3 (g) = g(fi + A 1 g 1 ' 2 ) ni +A 29 y 1 (133) 
gives the strong-coupling amplitude 

B 3 = -j- = 4.2835 , (134) 
A 2 

with an error e(B 3 ) = 48%. Certainly, such an accuracy is not sufficient. 

9.2 Iterated root approximants 

Using the iterated root approximants of Sec. 2, with setting a 6 = 0, we get the strong- 
coupling amplitudes 

B 2 = 8.713 , B 3 = 4.765 , B 4 = 3.2924 . (135) 

The error of the highest-order approximation is ^(-84) = 0.08%. This is the best accuracy 
among the considered approximants. 
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9.3 Inversion of strong-coupling series 



It is worth noting that in those cases when several terms of the strong-coupling expansion are 
known, it is possible to change the variable and to treat the old strong-coupling expansion 
as a weak-coupling one in terms of the new variable, and vice versa, the old weak-coupling 
expansion, as a strong-coupling one. For instance, the strong-coupling expansion for the 
ground-state energy of the Lieb-Liniger model can be written (see Ref. [41]) in the form 

, x 7T 2 / 4 12 32 80\ 

3 v g g 2 g 3 g A J 

Changing the variable to x — 1/ g makes the above series a weak-coupling expansion in terms 
of x. Conversely, Eq. (131) transforms, in terms of x, into a strong-coupling expansion [15]. 
However, possessing quite a number of terms of the strong-coupling limit with respect to the 
coupling parameter g is a rather rare case. Therefore, in the present paper, we focus on the 
problem of dealing with the direct coupling parameter, without changing the variable. 



10 Bose-Einstein condensation temperature 

An important problem is the dependence of the Bose-Einstein condensation temperature T c 
on the gas parameter 

7 = P 1/3 a s , (136) 

in which p is the average atomic density and a s , scattering length. For the asymptotically 
weak interaction, when 7 — > 0, the relative shift of the critical temperature, with respect to 
the condensation temperature of the ideal Bose gas, 



_ 2nh 2 
mkn 



P 



1 2 /3 



C(3/2) 



is linear in the gas parameter: 



AT 

V - c ^ (7 -> 0) , (137) 







where 



AT, = T r -T n 



The coefficient C\ has been calculated by a number of different methods (see review 
articles [42,43]). It can be represented as an asymptotic expansion 

Ci(g) ^ a x g + a 2 g 2 + a 3 g 3 + a 4 g 4 + a 5 g 5 (138) 

in powers of an effective coupling parameter [44-46], with the coefficients 

ai = 0.223286 , a 2 = -0.0661032 , a 3 = 0.026446 , 

a 4 = -0.0129177 , a 5 = 0.00729073 . 
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The sought value of c\ is the limit 

ci = lim a(g) = B . (139) 

The most accurate calculations of c\ have been accomplished by means of the Monte Carlo 
simulations in Refs. [47-51] and, by using different variants of the optimized perturbation 
theory [52], in Refs. [44-46,53-55]. Summarizing these results, we have: 

ci = 1.32 ±0.02 , [47,48] 

c 2 = 1.29 ±0.05 , [49,50] 

ci = 1.32 ±0.14 , [51] 

ci = 1.27 ±0.11 , [45,46] 

d = 1.30 ±0.03 . [54,55] 

The use of Pade approximants, for summing the asymptotic series (138), does not lead 
to accurate approximations, the best approximant having an error of —24%. 

10.1 Self-similar factor approximants 

The direct application of the self-similar factor approximants does not provide high accuracy. 
Thus, the factor approximant 

/*(«?) = a l9 (l + A l9 r(l + A 29 r , (140) 

under the condition 

1+711+712 = 0, (141) 

gives the amplitude 

B 3 = cuA^A? = 1.025 , (142) 
which approximates limit (139) with an error e(B 3 ) = —21%. 

10.2 Simple root approximants 

Simple root approximants also are not accurate enough. For example, the root approximant 

r* 2 (g) = a 1 g(l + A 1 g)- 1 (143) 

yields the amplitude 

B 2 = ?L = - ^ = 0.754 , (144) 
Ai a 2 

with an error e(B 2 ) = —42%. And the root approximant 

rl{g) = a l9 ((1 + A x g) ni + A 2 g 2 y l/2 (145) 

gives the amplitude 

B 3 = -^L= = 0.916 , (146) 
v A 2 

whose accuracy is only sightly better, having an error e(B 3 ) = —30%. 
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10.3 Iterated root approximants 

The iterated root approximants of Sec. 2 result in the amplitudes 

B 1 = 0.754 , B 2 = 1.383 , B 3 = 0.854 . (147) 

The last of these approximations has an error e(B 3 ) = —34%. The iterated root approximant 
of the fourth order does not exist, being complex valued. 

10.4 Corrected root approximants 

For the third-order correction function 

aig + a 2 g 2 + a 3 g 3 
" ^(J) (9^0), (148) 

the corresponding root approximant is 

C;(g) = l + b 2 g 2 (l + D ig )~ 2 . (149) 
The corrected root approximant becomes 

r 3 (9) = r* 2 (g)C*(g) . (150) 

This gives the amplitude 



^-|( 1+ !^ - 924 - < 151 » 



with an error £(-83) = —29%. 

Increasing the order of the correction function essentially improves the accuracy. Thus, 
for the fifth-order correction function 

a x g + a 2 g 2 + a 3 g 3 + a 4 g 4 + a 5 g 5 

CM " ' (152) 

whose root approximant is 

C*(g) = 1 + b 2 g 2 ((1 + Dig) 2 + D 2 g 2 )' 1 , (153) 
we have the corrected root approximant 

h(9) = r* 2 (g)C* 5 (g) . (154) 

The latter yields the amplitude 

5 5 = -^fl + -^-Ul.29, (155) 



a 2 \ Df + D 
which is in the frame of the numerically calculated values of c\. 
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10.5 Corrected-iterated root approximants 



The iterated root approximants, considered in subsection 10.3, are not accurate enough. But 
it is possible to correct them in the spirit of Sec. 2. Let us take, for example, the third-order 
iterated root approximant 

R 3 (g) = a l9 ((1 + A l9 f + A 2 g 2 y 1/2 . (156) 
Introduce the correction function 

a x g + a 2 g 2 + a 3 g 3 + a 4 g 4 + a 5 g 5 

cdg) " m ■ (157) 

whose self-similar approximant is 

C£(g) = 1 + 6 3 /(l + ^)- 3 . (158) 
The corrected-iterated root approximant 

Rl(g) = R 3 (g)C;(g) (159) 

gives the amplitude 

^ = ai(^ + ^)' 1/2 (l + ^) =1.31, (160) 

which is in perfect agreement with numerical calculations for c±. 

The fifth-order corrected root approximants of subsections 10.4 and 10.5 provide very 
accurate approximations for the coefficient c±. 



11 Conclusion 

In this paper, we have addressed the problem of extrapolating functions, which can be found 
only in the region of asymptotically weak coupling parameters, to the limit of infinitely strong 
coupling parameters. The main aim has been, assuming that the strong-coupling exponent 
is known, to find the strong-coupling amplitude. Several examples are considered, for which 
Pade approximants do not provide good accuracy or are not applicable at all. For these 
problems, we employ the self-similar approximation theory, comparing different variants of 
self-similar approximants. Two new classes of self-similar approximants are suggested, the 
corrected approximants and iterated approximants. 

Comparing the accuracy of different self-similar approximants for different problems, we 
find the following. For the Debye-Hiickel function (Sec. 3), the highest accuracy is provided 
by the iterated root approximants. For the quartic anharmonic oscillator (Sec. 4), the best 
are the self-similar factor approximants. For the polymer chain (Sec. 5), all considered 
variants of the self-similar approximants provide almost the same accuracy. For the oscillat- 
ing fluid string (Sec. 6), the best are the corrected root approximants. For the fluctuating 
fluid membrane (Sec. 7), the corrected root approximants are also the most accurate. For 
strongly interacting fermions (Sec. 8), the best are the corrected factor approximants. For 
the one-dimensional Bose system (Sec. 9), the iterated root approximants can give the best 
accuracy, though with a large dispersion between the different-order approximants. For the 
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Bose-Einstein condensation temperature (Sec. 10), the corrected root approximants give the 
best accuracy, yielding the results that are in very good agreement with the known numerical 
calculations. 

The considered examples illustrate that the self-similar approximants provide a versatile 
tool for extrapolation problems. 
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